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Abstract 

We use a local theory of photons purely as particles to model the 
single-photon experiment proposed by Tan, Walls, and Collett. Like 
Tan et al. we are able to derive a violation of Bell's inequalities for pho- 
ton counts coincidence measurements. Our local probabilistic theory 
does not use any specific quantum mechanical calculations. 
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1 Introduction 

The present paper is part of a research program which concerns a founda- 
tional analysis of phenomena usually described by quantum electrodynamics 
(QED) H Q |l0| . Our previous papers give a particle theory for diffrac- 
tion of light and the Casimir effect. The present paper is focused on another 
foundational topic. It remains to be seen how far the program we have 
undertaken can be carried. 
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A probabilistic theory of photons with well-defined trajectories is as- 
sumed. The wave properties come from the expectation density of the 
photons. The photons are also regarded as virtual, because they are not 
directly observable, including their annihilation of each other (see assump- 
tions bellow). What can be detected is the interaction with matter. The 
meaning of virtual used here is not the same as in QED. In summary, our 
basic assumptions are: 

• Photons are emitted by harmonically oscillating sources; 

• They have definite trajectories; 

• They have a probability of being scattered by matter; 

• Absorbers, like sources, are periodic; 

• Photons have positive and negative states (+- photons and — photons) 
which locally interfere, when being absorbed. 



by 



The expected density of ±-photons emitted at t in the interval dt is given 

s±(t) = ^(licos^t), (1) 

where to is the frequency of a harmonically oscillating source, A s is a constant 
determined by the source, and t is time. We used \ ± \ cos(ut) rather than 
cos(wt), to have a density that is nonnegative for all t and is between and 
1. If a photon is emitted at t' , < t' < t, then at time t the photon has 
traveled (with speed c) a distance r, where 

t-t' = - c (2) 

The conditional space-time expectation density of ±-photons for a spheri- 
cally symmetric source with given periodicity u is: 

h ± = 8^( 1 ±™ s "(t -"-)), ( 3 ) 

where A is a real constant. 

The scalar field defined in terms of the expectation density h±(t,r\u) is 
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where £q is a scalar physical constant. Using (|3|), (Q) may be rewritten for 
a spherically symmetric source as 



£ = £ °\4^ COSUJ ( t ~^)- (5) 



Applying the standard definition of average intensity, we get the expected 
result 

£%A 

8vrr 2 

Note that the standard bracket notation is used for time averaging, i.e., 
taking an expectation with respect to t. 

Since the absorber, or photodetector, behaves periodically with a fre- 
quency u, the probability px of absorbing a photon in detector X is given 
by 

Px = -(l + cos(a;t + V0), (?) 

where tp is an arbitrary phase that can be randomized. 

The expected number E t (X±) of each type of photon absorbed by de- 
tector X is the time-averaged product 

E t {X±) = {hZ(0)pxW), (8) 

where h±(9) is the expected density of photons (with a phase 9), and ip is a 
phase on detector X. Note that Et(X+), for example, is a random variable 
that is a function of 9 and tp. When we take expectation with respect to 
the distribution of 9 we use subscripts to make clear that the expectation 
is with respect to 9. The averaging is required because an absorption of an 
individual photon by an atom of a photodetector takes on average several 
orders of magnitude longer than the mean optical period of the photons, 
both theoretically and experimentally Q. 

As we previously assumed, during the process of absorption, photons 
with different states (positive and negative) annihilate each other. So, the 
expected number of photons to be detected in each detector X is: 

E t (X) = \E t (X+)-E t (X-)\. (9) 

We present here a violation of Bell's inequalities [|| § || with a local 
description of photons. 
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2 Experimental Configuration 

We are interested in the experimental setup proposed in |ll| and also dis- 
cussed in [12|| . The scheme uses two coherent sources ai(#i), with phase 9±, 
and 02(62)1 with phase 62, and a third source to be studied, u(8), with un- 
known phase. The experimental configuration has two homodyne detections, 
(D\, D2) being one and (-D3, -D4) the other, such that the measurements are 
sensitive to phase changes in u(6). The geometry of the setup is shown in 
Figure 1. In Figure 1 BS1, BS2 and BS3 are beam splitter mirrors that will 




Figure 1: Proposed experimental configuration. 

reflect 50% of the incident photons and let 50% of them pass. When pho- 
tons are reflected, the mirrors add a phase of tt/2 to the expected density, 
while no phase is added to the expected density when photons pass through 
BS1, BS2 or BS3. It is easy to devise a way to have the expected density 
of photons changed by a tt/2 phase by just delaying the photons that are 
reflected, and hence have interacted with the mirror, by a time T/4, where 
T is the period of the photon source. We will look for correlations between 
the pairs of photon detectors (D\,D2) and (D^^D^). 

The expected density of ±-photons, generated by the source u(0) is: 

fc£(0) = |(l±cos(wt + 0)). (10) 
The expected density coming from u(0) at each detector is: 

h^(e) = ^(i±co S Lt + 9 + ^\Y (11) 
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h^ 2 (e) = ^(i±cos(ut + e+Tr)), (12) 

8 

hg 3 {6) = - (l±cos(ut + e + ~)) , (13) 



(9) = (1 ± cos (ut + 9)). (14) 
8 

Note that we neglected factors of the form k • x coming from path contri- 
butions to the phase. We can do so considering the problem completely 
symmetric and remembering that only phase differences are relevant for the 
measurements we are proposing. 

In similar fashion, the expected density of ±-photons generated by the 
coherent sources oti, with phase 9{ and amplitude a/2, and ay, with phase 9j 
and amplitude a/2, in each detector, is given by the following expressions: 

h± 1 (0i) = j (licoefwt + fli + l)), (15) 

hg*{9 i ) = j(l±cos{ut + 9 i )), (16) 
h± z {9 j ) = |(l±cos(u;t + i )), (17) 

hi 4 (0j) = ~ (l ± cos (ut + % + |)), (18) 

where we again neglected path contributions to the phase and considered 
only the relevant phase at the detectors. 

We should point out that in equations (10) — (18) a and are split in 
half at each semi-mirror, because each time a photon reaches a mirror there 
is a probability of 1/2 that the photon passes through and a probability of 
1/2 that the photon is reflected by the mirror. 

The probability of absorption in each detector, consistent with equation 
([?]), is given most simply by the following equations. Some alternatives are 
formulated in equations (40) — ( p7|) at the end of this section. 

PDl = j(2 + cos (ut + 0i + + cos {ut + e + ^j, (19) 

C 

p D2 = — (2 + cos (ut + Oi + 7r) + cos (art + 9)) , (20) 

PA, = - ^2 + cos (ut + 9j) + cos (ut + 9 + V (21) 
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c 



PD 4 



2 + cos [tot + Oj + 



7T 



+ cos (ujt + 0)J 



(22) 



where C is a constant that corresponds to the efficiency of the detection 
process. 

The expected number of ± photons in each detector is given, according 
to equation (§), by the following expressions: 



Et(Df) = ( 






PDi) , 


(23) 


Et(Df) = < 




+ t&(0)) 


PD 2 ) > 


(24) 


Et(Dt) = < 




+ hg(0)) 


PD 3 ), 


(25) 


Et(Df) = ( 




+ h£*(0)] 


)PD 4 )- 


(26) 



Equations (23) — (26) use the fact that the expected number of photons at a 
detector is simply the sum of the number of photons from all sources. Also, 
in the equations above (F(t)} = y J$ F(t)dt, represents a time average of 
the random variable F(t), where t is time and (0,T) is a time interval such 
that ojT ^> 1. It is straightforward to obtain the expressions for the total 
expected number of photons in each detector from equations (||), @), fllq) — 
(0), ©— (H), and (H)— (HP, which we write as I k , for k = 1, . . . , 4, with 



Ik a function of 9 


and 9i or t 










h = 


Et(Dt) 


~ Et(Di) 


C 
" 16 


« + - 


+ (a + i/3) cos(# - 0;) 


, (27) 


h = 


Et(D+) 


-E t (D 2 ) 


C 
" 16 


a + - 


- (a + cos(# - e<) 


, (28) 


h = 


Et(D+) 


-Et(Dz) 


C 
~ 16 


a H — 
2 


- ( a + i/3)sin(0-%) 


. (29) 


h = 


E t (Dt) 


~E t {Dl) 


C 
~ 16 


° + 2 


+ (a + i/0) sin(6>-%) 


• (30) 



The expressions on the right hand side of (27) — (|30| ) are nonnegative, in- 
dependent of taking their absolute value, and so we subsequently drop the 
absolute values. 

We are interested in the correlation between the two pairs of detectors. 
First we need the variances 



Var e (Ii - h) = E e ((h - I 2 f) - (E e (h - I 2 )f, 



(31) 
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Var e (/ 3 - h) = E e ({h ~ hf) - (E e (h ~ h)f, (32) 
and covariance 

Cov ((/i-I 2 )(J 3 -/4)) = Ee((h-h){h-h))-Ee(h-h)Ee(h-h), (33) 

where Eq(I^) = ^- J Q 27r Ikd6, for k = 1, . . . , 4, is an expectation with respect 
to 9, with 9 uniformly distributed on [0, 2ir]. Thus 

Vareih ~ h) = ^C 2 ((3 + 2a) 2 , (34) 
512 

Var,(/ 4 - 7 3 ) = ^C 2 ((3 + 2a) 2 , (35) 
512 

and 

Cov„((Ii - I 2 )(h ~ h)) = ~^f 2 ^ + 2 «) 2 sin (^ " (36) 
The correlation is given by 

' ' J\<x,(h-h)\<n,(h-Uy 1 ' 

or 

Pe (h - I 2 , h ~h) = - sin(0i - d 3 ). (38) 

It is easy to show that (^) violates Bell's inequalities when four appropriate 
phases are chosen. 

An examination of the derivation of (|3^) shows that without serious 
change it holds simply for a classical field as (||). Details and discussion 
can be found in ||. In the case of both Varg and Covg, it is important to 
note that if we computed the correlation with respect to t rather than 9, we 
would get different results. It is easy to show, for example, that 

Var e E t (Df) ± Vsx t E e (Df). (39) 

In contrast, the order of 9 and t does not matter in analyzing the data of 
discrete photon counts, in Section 3. 

An attentive reader may object to our expression for the probability of 
detection, because we assume that the detector has the same probability to 
oscillate in phase with the noncoherent source as it has to oscillate with the 
coherent source, and that may bring some non-local characteristics to the 
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model. We can respond to this by examining the following probability for 
absorption. 

c r ( „ vr N 



PD 1 =j\l+COB\wt + 9i + -)), (40) 

C 

PD 2 = - (I + cos (tut + 6i)), (41) 

PD3 = |(l + cos(wt + ^)), (42) 



The probabilities above have no term on 9, but depend only on the phase of 
the coherent sources. This fact have the effect of wiping out all influences 
that the non-coherent source have on the detectors, and hence putting only 
local parameters, like 9i or 9j, depending on the detector, in the probability 
for detection. If we redo the computations for the correlation with the 
probabilities above, we end up with the same correlation function for a pair 
of homodyne detections. In fact, to point out the robustness of the result in 
face of the choice of probability for detection, we may examine the following 
set of probabilities. 

. acos(wt + 0i + f) + (3 cos (Lut + 9 + %)\ 

P*=o(l + ^ 2 -)> ( 44 ) 

. ,' a cos (cot + 9i + 7r) + (5 cos (cot + 9) \ 
^ = 2 I 1 + —p J ' (45) 

C A acos(ojt + 9j) + (3cos(ujt + 9 + f)\ 
PD3= 2l, 1 + aTT^ J' (46) 

CA acos(cjt + % + §) +/3cos(wt + e)\ 
TO *=2( 1 + " ^ J' <47) 

The above expressions would have a different physical interpretation from 
the previous two presented. Each phase is given a probability that is pro- 
portional to the amplitude of the source with the corresponding phase. The 
stronger the source, the more probable to find the detector with the same 
phase. It is again easy to show that if we use these probabilities we get the 
same correlations as before. 
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3 Photon Counts that Violate Bell's Inequalities 



In this section we are going to use the previous result to model discrete 
photon counts in such a way that they violate Bell's inequalities. For this, 
we define two new discrete random variables X = ±1 and Y = ±1. These 
random variables correspond to nearly simultaneous correlated counts at the 
detectors, and are defined in the following way. 

y _ \ +1 if a photon is detected at D\ . . 

1—1 if a photon is detected at D2 

Y _ f +1 if a photon is detected at D 3 , . 

J — 1 if a photon is detected at D4. 

To compute the expectation of X and Y we use the stationarity of the 
process and do the following. First, let us note that 

h-h = N x -P(X = l)-N x -P(X = -1), (50) 

where Nx is the expected total number of photons detected at D\ and D2 
and P(X = ±1) is the probability that the random variable X has values 
±1. The same relation holds for 

I 3 -I i = N Y P(Y = 1)-N Y -P(Y = -1). (51) 

To simplify we put as a symmetry condition that Nx = Ny = N, i.e., the 
expected number of photons hitting each homodyne detector is the same. 
But we know that 

I 1 + I 2 = N-P(X = l)+N-P(X = -l) = N, (52) 

and 

J 3 + J 4 = N ■ P(X = 1) + N ■ P(X = -1) = N. (53) 
Then we can conclude from equations (|27|) — (l30|) that 

E d (X\6) = ^— ^ = cos(# - 0i), (54) 
n + ±2 

Ed{Y\0) = I 4^ L = sm(e-e j ), (55) 

-<3 + -<4 
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where E d represents the expected value of the counting random variable. It 
is clear that if 9 is uniformly distributed we have at once: 

E{X) = E e {E d {X\9)) = 0, (56) 

E(Y) = E e {E d {X\e)) = 0. (57) 
We can now compute Cov(X,Y). Note that 

Cov(X,Y) = E(XY) - E(X)E(Y) 

= E e {E d {XY\9)) - E e {E d {X\9))E e {E d {Y\9)) (58) 

and so 

1 f 27T 

Cov(X,Y) = — / E d (XY\9)d9 
^ Jo 

~ I** E d {X\9)dO x -L r E d (Y\0)dO. (59) 

27T Jo Z7T Jo 

In order to compute the covariance, we also use the conditional independence 
of X and Y given 9, which is our locality condition: 

E d (XY\9) = E d (X\9)E d (Y\9), (60) 

because given 9, the expectation of X depends only on 9i, and of Y only on 
9j. Then, it is easy to see that 

p(X, Y) = Cov(X, Y) = - sin(0i -Oj). (61) 

The correlation equals the covariance, since X and Y are discrete ±1 random 
variables with zero mean, as shown in ( |56|) and (|57|), and so Var(X) = 
Var(y) = 1. It follows at once from ( |61~D that for a given set of 9{'s and 9j's 
Bell's inequalities are violated. 
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